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Abstract. I review here some recent result on the thermal conductivity 
of chains of oscillators whose hamiltonian dynamics is perturbed by a 
noise conserving energy and momentum. 



1. Introduction 

Let us consider a 1-dimensional chain of oscillators indexed by x G Z, 
whose formal Hamiltonian is given by 

^2 



(1.1) 



where qx indicate the displacement of the atom x from its equilibrium po- 
sition, and Px its momentum (we fix for the moment all masses equal to 
1). The potential V and W are some smooth positive function growing at 
infinity fast enough. The W potential is often called pinning. 

We want to understand the macroscopic properties of energy transport 
for the corresponding hamiltonian dynamics 

qx = Px, Px = -dq^T-t. (1.2) 

When we say "macroscopic energy transport" we are meaning a certain non- 
equilibrium evolution that we want to observe in a space-time macroscopic 
scale that should be specified. This space-time scale can be typical of the 
model and the same model can have distinct macroscopic scalings under 
which it will behave very differently. For example in the case that W = Q 
(unpinned system) , total momentum is also conserved and hyperbolic scaling 
(in which space and time are scaled in the same way) is a natural one. In the 
hyperbolic scaling energy is carried around by the momentum of the particles 
and macroscopic evolution equation is given by the Eulcr equations. 

Another important scaling is the diffusive scaling (time scaled as square 
of space) where transport of energy happens by diffusion. The macroscopic 
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equation is usually given by heat equation. For example, when the pin- 
ning potential W is present, total momentum is not conserved and nothing 
moves in the hyperbolic scale. Energy will move on the diffusive space- 
time scale, and macroscopically its evolution should be governed by heat 
equation. More precisely, defining the empirical distribution of the energy 

£%G, t)=eY^ G{£x)£^{e-^t) (1.3) 

X 

2 

where Sx = ^ + V{qx+i — Qx) + W{qx) is the energy of particle x, and G 
is a smooth test function on M with compact support. We expect that 

lim£%G,t) = I G{y)u{y,t) dy (1.4) 

in some statistical sense, with u{y, t) solution of the (non-linear) heat equa- 
tion 

dtu = dx{K{u)dxu) . (1.5) 

The function k{u) is called thermal conductivity and can be expressed in 
terms of the dynamics in equilibrium: 

"^"^ = Bl. 2^ ^ [(f.(t)f (0))^ - u^] (1.6) 

X 

where the temperature T = T{u) correspond thermodynamically to the 
average energy u, < ■ >t is the expectation of the dynamics in equiHbrium 
at temperature T. 

Proving (1.4)- (1.5) from an hamiltonian microscopic dynamics is one of 
the major challenge in non-equilibrium statistical mechanics [7]. It is not 
clear under which conditions on the interaction and initial conditions this re- 
sult could be valid. Even the proof of the existence of the limit (1.6) defining 
the thermal conductivity is completely open for any deterministic system. 
What is clear is that (1.4)-(1.5) are not always valid. For example if V and 
W are quadratic (harmonic chains) , then the energy corresponding to each 
Fourier mode is conserved and carried ballistically without any interaction 
with the other modes. It results that thermal conductivity is infinite in this 
case ( [16], and for a macroscopic equation in the hyperbolic scaling see [8], 
as explained in section 4). 

In nonlinear unpinned cases {W = 0) one expects, generically in dimen- 
sion 1, that K = -)-oo, and correspondingly a superdifTusion of energy. This 
fact seems confirmed by all numerical simulation of molecular dynamics [11], 
even though there is not a general agreement about the order of this supcrdif- 
fusivity. Most interesting would be to understand what kind of stochastic 
process would govern this superdiffusion. 

Adding a stochastic perturbation to the hamiltonian dynamics certainly 
helps to obtain some mathematical result in these problems. Of course 
adding noise to the microscopic dynamics may change the macroscopic be- 
havior. The ideal is to add noise terms that change as little as possible the 



ENERGY DIFFUSION 



3 



macroscopic behaviour, at least qualitatively. For example it is important 
that this noise conserves energy, and eventually momentum, since these are 
the quantities we expect being conserved in the infinite system. 



2. Conservative Stochastic Dynamics 

We consider the Hamiltonian dynamics weakly perturbed by a stochastic 
noise acting only on momenta and locally preserving momentum and kinetic 
energy. The generator of the dynamics is 

L = A + jS (2.1) 

with 7 > 0, where A is the usual Hamiltonian vector field 



(2.2) 

while S is the generator of the stochastic perturbation. The operator S 
acts only on the momenta {py} and generates a diffusion on the surface of 
constant kinetic energy and constant momentum. S is defined as 

^=^E(^-)'' (2-3) 

where 

Yz = {Pz - Pz+i)dp^-i + {Pz+i - Pz-i)dp^ + {pz-i - Pz)dp^+i 

which is a vector field tangent to the surface of constant kinetic energy and 
of constant momentum for three neighbouring particles. As a consequence 
energy and momentum are locally conserved which, of course, implies also 
the conservation of total momentum and total energy of the system, i.e. 
formally 

sY,Px = o, sn = o. 

Since also the hamiltonian dynamics conserves energy, we also have LTi = 0. 
Furthermore in the unpinned case {W = 0), L X^^-Pa; = 0. 

The evolution of {p{t),q{t)} is given by the following stochastic differential 
equations 

dqx = Px dt, 

dpx = - dxHdt + - A(4p^ + px-i + Px+i)dt 

6 (2.4) 



+ V3 XI {Yx+kPx)dW:c+k{t)- 
* fc=-l,0,l 



Here {wy{t)}y^z are independent standard Wiener processes and A is the 
discrete laplacian on Z: A/(z) = f{z + 1) + f{z — 1) — 2/(z) . 
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In the unpinned 1-dimensional case, the equihbrium measures are par- 
ticularly simple. In fact the right coordinates are = Qx+i — Qx, and the 
family of product measures 

^.,,^Adp,dr) = ll—^ (2.5) 

are stationary for the dynamics. The three parameters A, tt, P correspond to 

the 3 conserved quantities of the dynamics (energy, momentum and '^^Tx, 
the stretch of the chain), while Z{X,P) is the normalization constant. It can 
be proven that these are the only translation invariant stationary measures 
of the dynamics ( [13], we call this property "ergodicity of the infinite dy- 
namics" ) . In more dimensions the unpinned case is much more complex and 
even the definition of these equilibrium measures are problematic (cf. [9]). 

In the unpinned one-dimensional case, since momentum is conserved, 
there is a non trivial macroscopic evolution on the conserved quantities 
in the hyperbolic scaling (space and time scaled in the same way). Let us 
define the energy of particle x as 

= f + ^ {V{rx-i) + V{rx)) 
Locally the conservation of energy can be written as 

LSx = jx-i,x ~ ix,x+i (2-6) 

where the instantanueous current of energy j^x+i ^he sum of the current 
due to the hamiltonian mechanism plus the current due to the stochastic 
term of the dynamics: 

3x,x+\ Jx,x+l ' 3x,x+l 

sa 1 , es (2.7) 

with = {pl^i + Apl + pI_^ + Px+iPx-i - '^Px+iPx - '^PxPx-i)- Similarly 
conservation of momentum reads as 

Lpx = jx-l,x ~ 3x,x+l^ 

fx,x+i = y'irx) + ^ViApx + px-i + Px) 
while mass conservation is simply given by 

Ltx = Px-^i - Px- (2.9) 

Let G{y) a test function continuous with compact support. We expect 
that 

eY,G{ex)lpx{e-H)] ^^^'"^ G{y) I p{t,y)] dy (2.10) 
\Sx{e-H) Vit,y) 
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where r(t, y),p{t, y), e(t, y) are given by the solution of the Euler hyperbolic 
system of equations 

dtt = dyp 

dtP = dyPit,t-py2) (2.11) 

% = 5^(pP(r,e-pV2)) 

Here P{r, u) is the thermodynamic pressure, which is related to the thermo- 
dynamic entropy S{r,u) by the relation 

and S is defined from Z(X, (3) with a Legendre transform: 

S(r, u) = sup iXr - Pu - log (z{X, P)y/p/2Tr] } (2.13) 

Pressure P{r,u) is also give by the expectation of V'{rx) with respect to 
/^A,7r,/3) foi' the corresponding values of the parameters A and (3. 

The hydrodynamic limit (2.10) can be proven rigorously in the smooth 
regime of equation (2.11), by using the relative entropy method ( [14], see 
also [1] and [13] for the application to this specific model). Note that (2.10) 
does not depend on the strength of the microscopic noise 7. Noise here 
is used only to prove some ergodic properties for the dynamics, necessary 
to obtain the result. In fact without noise this limit may not be true, as 
for example in the linear case (V quadratic, see below). Note also that for 
smooth solutions the macroscopic evolution (2.11) is locally isoentropic, i.e. 

|s(t(t,y),e(t,y)-P^)=0 (2.14) 

A challenging open problem is to extend this result to solutions that present 
shocks, where the above derivative is (presumably) strictly positive. 

3. Diffusive evolution: Green-Kubo Formula 

In the pinned model {W > 0) momentum is not conserved, energy is the 
only relevant conserved quantity for the infinite system and the equilibrium 
measure arc the Gibbs measures at given temperature, corresponding to the 
hamiltonian 7i. These probability measures are defined by the usual DLR 
equations. Consequently the hyperbolic scaling is trivial (nothing moves at 
that time scale). In order to see energy moving at a macroscopic scale, one 
has to look at larger time scale. The natural scaling is the diffusive one: we 
expect, for a given test function G as above, 

e ^ G{ex)£Ae-'t) / G(y)T(i, y) dy (3.1) 

where 'Z{t, y) is the solution of the (non-linear) heat equation 

dtl = {K{%)d^V) , (3.2) 
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where n(T) is the thermal conductivity at temperature T. This is given by 
the Green-Kubo formula 

+00 



2x(T) 



poo J ^ 



(3.3) 



where (|jf '"+i(i)io'i'(^)^ denote the expectation with respect to the dy- 
namics in equilibrium at temperature T = j3^^. The explicit 'jT'^ /2x{T) 
term is the contribution of the stochastic part of the dynamics. Formula 
(3.3) can be obtained from (1.6) using the conservation of energy (see [4] 
for a proof). We believe that such statement is always true for non-linear 
pinned dynamics, also in the deterministic case (7 = 0). But even for 7 > 0, 
this is hard to prove and still an open problem. Even the convergence of 
the integrals defining (3.3) is not known if non-linearities are present (some 
bounds are proven in [4]). The only case with non-linear interaction in 
which (3.3) is proven convergent is when the noise is generated by Langevin 
heat bath attached at each particle of the system [6], a non conservative 
stochastic perturbation. 

In unpinned systems the situation is more complex because of the momen- 
tum conservation. To avoid complicate re-centering along characteristics of 
(2.11), we can consider initial random configurations with momentum of 
(locally) zero average and constant density profile, only gradients of tem- 
perature admitted. Then conductivity k is also a function of the density of 
particles and in formula (3.3) the expectation should be taken with respect 
to the equilibrium dynamics with tt = (i.e. starting with configurations 
distributed by Ma,o,/3 defined by (2.5)). Numerical evidence shows that in 
this one-dimensional case k = +00 ( [11], also for 7 > 0, as long as mo- 
mentum is conserved [2]). In the physics literature there is a long discussion 
about the nature and the order or this superdiffusion. From dimension 3 on, 
it is expected that formulas corresponding to (3.3) give a finite diffusivity. 

Rigorous results can be proven for the harmonic case with 7 > [3,4]. It 
turns out that k is finite if system is pinned or in dimension d > 3, while 
the 1 and 2 dimensional unpinned cases are super diffusive. 



4. Kinetic Limits: Phonon Boltzmann Equation 

The harmonic case is enough simple to obtain some non-equilibrium re- 
sults. In [5] we consider the hyperbolic scaling in a weak noise limit. Noise is 
rescaled by multiplying its strength 7 by e. This way the effect of the noise 
per particle remains finite in the macroscopic scale (that motivates the term 
kinefAc in defining this limit). This is in the same spirit as the model of hard 
sphere with random collision considered in [15]. The right quantity to look 
here is the Wigner distribution of the energy, formally defined as 

W'{y, k,t) = - / e'^'^yi/' {i/jik - r]/2, t/e)*i{;{k + r]/2, t/e)) dr] (4.1) 

2 7-1/2 
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where 

= ^ i.oj{k)q{k, t) + i'Pik, t)) , 

here q{k,t),p{k,t) are the Fourier transform of the qy{t),py{t), and uj{k) is 
the dispersion relation of the lattice, which in this one-dimensional nearest 
neighbour case is given by uj{k) = c|sin(7rA;)| {acoustic dispersion). The 
result in [5] states that k, t) converges, as a distribution on M x [0, 1], 

to the solution of the linear transport equation 

dtWiy,k,t)+'^dyWiy,k,t) = ^ J Cik,k'){W{y,k',t)-Wiy,k,t)) dk . 

(4.2) 

In the deterministic case (7 = 0) this result was obtained by Dobrushin 
et al. in [8], see also [12]. The collision kernel C{k,k') is positive and 
symmetric. It is computable explicitly (cf. [5]), but the important point is 
that C{k,k') ~ k'^ for small k. This is a consequence of the conservation 
of momentum: long waves scatter very rarely. Because C{k, k') is positive, 
(4.2) has a simple probabilistic interpretation: W is the density at time t 
of the energy of particles (phonons) of mode k. A phonon of mode k moves 
with velocity uj'{k) and after an exponentially distributed random time of 
intensity jC{k, k') changes its mode to k' . Defining a Markov jump process 
K{t) in [0, 1] with jumping rate 7C, the position of the phonons is given by 
Y{t) = j',J{K{s))ds. 

Thermal conductivity can be computed from (4.2) (cf. [5]) and the results 
are in accord with the direct calculations done in [4]. 

5. Levy's Superdiffusion of Energy 

As we mention in the previous section, phonons of small k scatter rarely, 
but their velocity uj'{k), in the unpinned case, are still of order 1 as ^ 0. 
This induces a superdifFusive behavior of these phonons. In [10], as applica- 
tion of new limit theorems for functionals of Markov chains and processes, 
we prove that, for a = 3/2, 

6y(6-"t) c{t) (5.1) 

law 

where C{t) is a Levy a-stable process, i.e. a stochastic process with inde- 
pendent stationary increments and with C{1) distributed by a ct-stablc law. 
In terms of the solution W{y, k, t) of the equation (4.2) this result implies 
the convergence 

lini J \W{e-''t,e-^y,k) - u{t,y)f dk = (5.2) 

where u{t, y) is the solution of the fractional heat equation 

dtu = -c{-Ay)'^/\ . (5.3) 
where c is a positive constant. 
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In the pinned case all the above result are still valid, but since the velocity 
of the phonons u>'{k) ~ k for small k, we have a regular diffusive behavior, 
and a = 2. 

It would be very interesting to understand how these results extends to 
the anharmonic cases. Equation (4.2) will be substituted by the non-linear 
phonon Boltzmann equation [17]. In the unpinned one-dimensional case 
thhis equation still will produce a superdiffusion. Is it again of Levy type, 
or will have some non-markovian terms? 
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